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Abstract. We study in random-phase approximation the newly discovered supersolid phase of 4 He and 
present in detail its finite temperature properties. 4 He is described within a hard-core quantum lattice gas 
model, with nearest and next-nearest neighbour interactions taken into account. We rigorously calculate 
all pair correlation functions in a cumulant decoupling scheme. Our results support the importance of the 
vacancies in the supersolid phase. We show that in a supersolid the net vacancy density remains constant 
as function of temperature, contrary to the thermal activation theory. We also analyzed in detail the 
thermodynamic properties of a supersolid, calculated the jump in the specific heat which compares well to 
the recent experiments. 
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1 Introduction 

One of the biggest accomplishments of theoretical con- 
densed matter physics is the ability to classify various 
phases and phase transitions by their mathematical or- 
der. These mathematical orders are usually expressed by 
order parameters reflecting certain limiting behaviours of 
two particle correlation functions. This concept makes it 
possible to describe the properties of physically very dif- 
ferent systems in a common language and establish a link 
between them. 

Without this concept, the counter-intuitive idea of a su- 
persolid, i.e. a solid that also exhibits a superflow, would 
not have been conceivable. In this language supersolidity, 
firstly proposed by Andreev and Lifshitz[l in 1969 and 
picked up by Leggett and Chester in 1970 [3l[2] is the clas- 
sification of systems that simultaneously exhibit diagonal 
and off-diagonal long range order. Yet the idea of a super- 
solid was reluctantly received by the scientific community, 
because early experiments failed to detect any such effect 
in Hclium-4. Apart from John Goodkind 4 who had pre- 
viously seen suspicious signals in the ultrasound signals of 
solid helium, physicists were surprised when in 2004 Kim 
and Chan [5|[6] announced the discovery of supersolid he- 
lium. Although equipped with a head start of more than 
30 years the theoretical understanding of the supersolid 
state lags vastly behind, not least because supersolidity 
as we observe it today is nothing like what the pioneers 
in early 1970's had anticipated. 

It is now evident that defectons and impurities play a cru- 
cial role in the formation of the supersolid state. However, 
the data and the results of various experiments draw a 



rather complex picture of the supersolid phase. The nor- 
mal solid to supersolid transition is not of the usual first or 
second order type transition but bears remarkable resem- 
blance to the Kosterlitz Thouless transition well-known in 
two dimensional systems. Furthermore annealing experi- 
ments show that 3 He impurities play a significant role but 
the data remain all but conclusive, as the measured su- 
perfluid density varies by order of magnitudes among the 
different groups. Recently a change in the shear modulus 
of solid helium was found and the measured signal almost 
identically mimics the superfluid density measured by Kim 
and Chan [5H5]. Some popular theories give plausible ex- 
planations of some aspects of the matter. That, such as 
the vortex liquid theory suggested by Phil Anderson [9], 
is in good agreement with properties of the phase tran- 
sition; theories based on defection networks are capable 
of describing the change in the shear modulus. However, 
we feel that a satisfying and comprehensive theory is still 
missing. Part of the problem stems from the complexity of 
the models. Many models are only at sufficient accuracy 
solvable with numerical Monte Carlo methods. The re- 
sults, doubtlessly useful, are seldom intuitive and the lack 
of analytical results do not meet our notion of the under- 
standing of a phenomenon. Other approaches on the other 
hand are limited in their analytical significance and rather 
represent a phenomenological ansatz. 
In this work we attempt to fill a part of this gap in the the- 
oretical understanding of the supersolid phase. We present 
a theory of supersolidity in a Quantum-Lattice Gas model 
(QLG) beyond simple mean-field approaches. Following 
the approach of K.S. Liu and M.E. Fisher [T2"] we map the 
QLG model to the anisotropic Heisenberg model (aHM). 
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The method of Green's functions proves to be very suc- 
cessful in the description of ferromagnetic and antiferro- 
magnetic phases and we use this method to investigate the 
supersolid phase which corresponds to a canted antiferro- 
magnetic phase. Applying the Random-phase Approxima- 
tion (RPA), we broke down using cumulant decoupling, 
the higher order Green's functions to obtain a closed set 
of solvable equations. This method gives a fully quantum 
mechanical and analytical solution of a supersolid phase. 
We will see that quantum fluctuations have a significant 
effect on the net vacancy density of the supersolid and we 
will also see that the superfluid state is unstable against 
zero-energy quasi-particle excitations. We also derive im- 
portant thermodynamic properties of this model and de- 
rive formulas for interesting properties such as the power 
law exponents and the jump of the specific heat across the 
normal solid to supersolid line. 

The paper is organized as follows: In Section [2] we intro- 
duce the generic Hamiltonian of a bosonic many body sys- 
tem and discretize it to a quantum lattice gas model. This 
model is equivalent to the anisotropic Heisenberg model 
in an external field and we will identify the corresponding 
phases in Section [3] In Section |4] we re-derive the classical 
mean-field solution and discuss briefly their significance 
before we in Section [5] recapitulate the Green's functions 
for the anisotropic Heisenberg model in the random-phase 
approximation at zero temperature. The ground state of 
the system is obtained by solving the corresponding self- 
consistency equation. In the following two sections we de- 
rive basic thermodynamic properties as well as ordinary 
differential equations to calculate the first and second or- 
der phase transition lines. In Section [8] we analyse the 
quasi-particle energy spectrum and in the last section we 
calculate phase diagrams for various parameter sets and 
analyse the properties of the supersolid phase and the cor- 
responding transitions. 



2 Generic Hamiltonian and Anisotropic 
Heisenberg Model 

If we neglect possible 3 He impurities supersolid Helium-4 
is a purely bosonic system whose dynamics and structure 
are governed by a generic bosonic Hamiltonian: 

H= [ d 3 x^(x)(——V 2 + fiU(x) 
J 2m 

+~ J d 3 xd 3 x'ip^(x)^(x')V(x-x')t/j(x)tlj(x') 

(1) 

Here tp\x) are the particle creation operators and ip(x) 
destruction operators and obey the usual bosonic commu- 
tation relations. Hamiltonians such as in Eq. fT]) are not, 
due to the vast size of the many body Hilbert-space, di- 
agonalisable even for simple potentials V(x). The accom- 
plishment of any successful theory is to find an approxi- 
mation that sustains the crucial mechanism while reduc- 
ing the mathematical complexity to a traceable level. Here 



we follow the method of Matsubara and Matsuda [H] who 
successfully introduced the Quantum Lattice Gas (QLG) 
model to describe superfluid Helium. We believe that the 
quantum lattice gas model is particularly useful for super- 
solids as the spatial discretization of this model serves as a 
natural frame for the crystal lattice of (super)-solid helium 
and a bipartite lattice elegantly simplifies the problem of 
breaking translational invariance symmetry for states that 
exhibit diagonal long range order. 

According to Matsubara and Tsuneto [13] the generic Hamil- 
tonian Eq. Jl} in the discrete lattice model reads: 

H = /i^rij + u ij( a i ~ a ]){ a i - a j) + VijUiUj 

i ij ij 

(2) 

Here Uij are solely finite for nearest neighbor and next 
nearest neighbor hopping and otherwise zero. The val- 
ues of u nn and u nnn are such that the kinetic energy is 
isotropic up to the 4th order. As atoms do not penetrate 
each other there can exist only one atom at a time on a lat- 
tice site and consequently <v and a are the creation and 
annihilation operators of a hard core boson commuting 
on different lattice sites, but obey the anti-commutator 
relations on identical sites. Equation Q is the Hubbard 
model in 3 dimensions for hard core bosons. Neither fully 
bosonic nor fermionic, the lack of Wicks theorem inhibits 
the application of pertubative field theory though hard 
core bosonic systems are algebraic identical to spin sys- 
tems. A simple transformation, where the bosonic opera- 
tors are substituted by spin-1/2 operators generating the 
SU(2) Lie algebra, maps the present QLG model onto the 
anisotropic Heisenberg model: 

h = h* ]T si + ]T 4s?s; + ]T 4(s?s* + sfs^s) 

i ij ij 

Here the correspondence between the QLG parameters 
and the spin coupling constants is given by: 

4 = ^ 

3 3 

3 Phases 

The self-consistency equations in the random-phase ap- 
proximation as we will derive in a later chapter are very 
lengthy and therefore it is our primary goal to keep the 
present model as simple as possible while still being able to 
describe the crucial physics. For this reason we will define 
the anisotropic Heisenberg model on a bipartite BCC lat- 
tice which consists of two interpenetrating SC sub-lattices 
as Figure [T] shows. In this lattice geometry the perfectly 
solid phase is composed of a fully occupied (on-site) sub- 
lattice A while sublattice B refers to the empty interstitial 
and is consequently vacant. As there is no spatial density 
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Fig. 1. The BCC lattice consists of two interpenetrating SC 
sub-lattices. In the perfect solid phase one sub-lattice (i.e. sub- 
lattice A) serves as on-site centers and is occupied while sub- 
lattice B represents the empty interstitial and is vacant (For 
simplicity, we only present the two dimensional case). 
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Fig. 2. The phase diagram for J? = 1.498, Jj = 0.562, 
jj 1 = -3.899 and jj = -1.782 as calculated by MF. 



Table 1. Possible magnetic phases and the corresponding 
phases of the Hubbard model. All Phases are defined by their 
long range order. The columns, from left to right, are the spin 
configurations, magnetic phases, ODLRO, DLRO and corre- 
sponding 4 He phases, respectively. 



TT 


FE 


No 


No 


Normal Liquid 


// 


CFE 


Yes 


No 


Superfluid 


// 


CAF 


Yes 


Yes 


Supersolid 


Tl 


AF 


No 


Yes 


Normal Solid 



variation in the liquid phases both sublattices are equally 
occupied, and the mean occupation number simply cor- 
responds to the particle density. We are aware that the 
above choice of SC sublattices does not reflect the true 4 He 
crystal structure which is hexagonal closed-packed (HCP). 
Nevertheless we believe that crucial physical properties 
such as phase transition and critical constants do not de- 
pend on the specific geometry as long as no other effects 
such as frustration are evoked. 

Table |3] charts the various magnetic phases and identifies 
the corresponding phases of the 4 He system. According 
to their spin configurations we identify the four magnetic 
phases to be of ferromagnetic (FE), canted ferromagnetic 
(CFE), canted anti- ferromagnetic (CAF) and antiferro- 
magntic (AF) orders. The off-diagonal long range order 
parameter m\ and the diagonal long range order mi are: 

mi = (S* A ) + (S%) 

m 2 = (S\) - (S" B ) (5) 

which readily identify the corresponding phases of the He- 
lium system. 



4 Mean-Field Solution 

In our previous work 10J we have already re-derived the 
classical mean-field solution of the anisotropic Heisenberg 
model at zero temperature. As this model provides an easy 
and intuitive access to the model we extend the approx- 
imation to finite temperatures as was done by K.S Liu 
and M.E. Fisher [12] and briefly discuss its solution and 
phase diagram. The mean-field Hamiltonian is obtained by 
substituting the spin-i operators with their expectation 
values: 

Hmf = -h'({S' A ) + {S' B )) 
-24(S A )(S B ) 4({S1)(S%) + (S' B )(S%)) 
-2J?(S%)(S%) - JJ((S A )(S A ) + (S%)(S%)) 

(6) 

Here j| = -q\J^ AtjeB , j\ = -q2jf eA j eA , J% = 
-q\JJ eA ,j^g and jj = -qiJj £A]£A where qi = 6 and 
<72 = 8 are the number of nearest and next nearest neigh- 
bors. The mean value of S y drops out as the randomly 
broken symmetry S x <-> S y (ODLRO) allows for (S y ) = 0. 
The standard method of deriving the corresponding self- 
consistency equations is to minimize the Helmholtz's Free 
energy F = H — TS. The entropy S is given by the pseudo 
spin entropy of the system: 

S = ~\[{\ + S A ) In(i + S A ) + (i - S A ) hi(\ - Sa) 

+ {\ + S B ) ln(i + S B ) + {\- S B ) ln(~ - S B )}(7) 

where S A = y/{S ZA )* + (S XA )* and S B = V \S ZB )^ + (S XB )* 
In the canted anti-ferromagnetic and the canted ferro- 
magnetic states there are 4 self-consistency equations in 
the ferromagnetic and anti-ferromagnetic phases where 
{S x ) = {S y ) = they are reduced by two. These equations 
are readily obtained by differentiating the free energy with 
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respect to (S z ) and (S x ) respectively. The resulting equa- 
tions can be rearranged to yield: 

V(Sza) 2 + (s XA ) 2 = t .^i^A 
v / (^FT(^=^f^(8) 



h* + 2{S ZA )(4 - Jj) + 2(S ZB )4 = 2Jl 



(S XB ) 
(S XA ) 



h* + 2(S ZB )(4 Jj) + 2(S ZA )4 = 2J T ^4-(S ZB )(9) 



where 



u lA = [(2J^(S B )+2jJ(S x A )f + 
(24(S z B )+24(S%) + h z ) 2 }i 

u B = [(2J^(S A )+2jJ(S B )) 2 + 
{2j\\{S B )+2jW{S z A ) + h z ) 2 ]? 



(10) 



The two equations (Eq. ([5])) are dismissed in the ferromag- 
netic and the anti-ferromagnetic phases as the transversal 
fields {S A ) and (S B ) become zero. At zero temperature 
T = the free energy and (H) coincide. This allows us to 
deduce the phases at absolute zero in the limiting cases 
(limits of h z ) from Equation ©. In the limit of h z — > oo 
the Hamiltonian reduces to an effective one particle model: 



H = h z ((S ZA ) + (S ZB )) 



(11) 



Consequently the system will be in the energetically favor- 
able ferromagnetic phase. In the opposite limit, h z — > 
and with antiferromagnetic coupling j| < the term of 
nearest neighbor interaction 



H — J\ (S ZA )(S ZB ) 



(12) 



is the only term that significantly lowers the energy. There- 
fore the ground-state is the anti-ferromagnetic state. At 
T=0 and a suitable choice of parameters the canted fer- 
romagnetic and the canted anti-ferromagnetic phases are 
realized in-between these limits as seen in Figure [5] How- 
ever with increasing temperature the regions of canted 
ferromagnetic and canted anti-ferromagnetic phases de- 
plete and at sufficiently high temperatures only the anti- 
ferromagnetic and ferromagnetic phases survive. As we 
can see from the phase diagram in Figure [2] most tran- 
sitions are of second order. Only for parameter regimes 
where the CAF does not appear the resulting CFE-AF 
is of first order. The boundary lines are defined by ordi- 
nary differential equations and generally need to be calcu- 
lated numerically. Nonetheless there exist analytical ex- 
pressions for the locations of the transitions at absolute 
zero as derived by Matsuda and Tsuneto [13]. The ferro- 
magnetic to canted ferromagnetic transition point is de- 
termined by Equations §§§ if we set (S ZA ) = (S ZB ) = \ 
and (S XA ) = (S XB ) = 0: 



l FE — CFE 



(13) 



Equally the canted anti-ferromagnetic to anti-ferromagnetic 
transition is defined by (S ZA ) — —(S ZB ) — \ and (S XA ) = 
(S XB ) = 0: 



L C AF—AF 



fi- 



4 + 4-jJ) 2 -(jI) 



(14) 



The canted ferromagnetic and the canted anti-ferromagnetic 
phases coexist where the order parameter of the DLRO 
,mi = (S ZA ) — (S ZB ) approaches zero. We replace (S ZA ) 
and {S Zb ) in Equation (J9j) with mi andm2 = (S Za ) + (S Zb ) 
and retain only linear terms of mi. Subtracting and adding 
both equations respectively yields: 



2mi(4' 



2m 2 {4 - 
J 2 ~ 4 , 



Jj 



4 



„ tT „ Am 2 
-2J 1 T ml 



2J{ m 2 
1 



4m,2 



1 



(15) 



We used that ^{S ZA ) 2 + (S XA ) 2 = \ at T=0. The solu- 
tion of these two equations determine the corresponding 
transition point which is given by: 



CFE—CAF 



J 1 + J 2 + J J 4 



(-4 + 4 



(16) 



5 Green's Functions 



Although the classical mean-field theory is quite insight- 
ful and gives an accurate description of the variously or- 
dered phases its fails to take quantum fluctuations and 
quasi-particle excitations into account. Hence, in order to 
overcome these shortcomings and to derive a fully quan- 
tum mechanical approximation we solve the anisotropic 
Heisenberg model in the random-phase approximation which 
is based on the Green's function technique. At finite tem- 
perature the retarded and advanced Tyablikov P~5l[T5] com- 
mutator Green's function defined in real time are: 



(t) = -ie(t){\[snt),s»] 
(t) = ie(-t)(\[s?(t) > s»] 



(17) 



The average () involves the usual quantum mechanical as 
well as thermal averages. The most successful technique of 
solving many body Green's function involves the method 
of equation of motion which is given by: 



idtG. 



IjAdv 



(t) 



= 5(t)([St,S?])-i6(t)([[S?,H] 
6{t){[S*,S?})+ie(-t)([[S!,H] 



(18) 



The commutator [Sf , H] can be eliminated by using the 
Heisenberg equation of motion giving rise to higher, third 
order Green's functions on the RHS. In order to obtain a 
closed set of equations we apply the cumulant decoupling 
procedure which splits up the third order differential equa- 
tion into products of single operator expectation values 
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and two spin Green's functions. The cumulant decoupling 
p"7] is based on the assumption that the last term of the 
following equality is negligible: 

(ABC) = 
(A)(BC) + (B)(Ad) 
+ (C)(AB)- 2(A) (B)(6) 
+ ((A-(A))(B-(B))(C-(C))) (19) 

The set of differential equations is not explicitly dependent 
of the temperature, i.e. temperature dependence solely 
comes with thermal averaging of the single operator expec- 
tation values. Therefore the solution is formally identical 
to the zero temperature solution and the detailed deriva- 
tion of the Green's function in their full form can be found 
in Ref. [HI- 

The averages of the spin operator, appearing in the cu- 
mulant decoupling scheme, have to be determined self- 
consistcntly. Two self-consistent equations can be derived 
from correlation functions corresponding to the Green's 
functions. The self-consistency equations of the canted 
ferromagnetic (superfluid) and canted antiferromagnetic 
(supersolid) phases are structurally different from the fer- 
romagnetic (normal solid) and antiferromagnetic (normal 
fluid) phases, as the off-diagonal long range order increases 
the number of degrees of freedom by two and hence two 
additional conditions, resulting from the analytical prop- 
erties of the commutator Green's functions, apply. Thus, 
the self-consistency equations of the canted phases, can 
be written as a function of the temperature, the external 
magnetic field and the spins in x-direction: 



A ((S A ),(S B ),h*,T) = 



F 



h z ,T) = 



(20) 



similar the self-consistency equations for the ferromag- 
netic and anti-ferromagnetic phases: 



F2 c ((S A ),(S B ),h*,T)=Q 
F% c ((S A ),(S* B ),h z ,T)=0 



(21) 



These equations involve a 3 dimensional integral over the 
first Brillouin zone. As explained in the work |10j on the 
zero temperature formalism this integral can be reduced 
to a two dimensional integral by introducing a general- 
ized density of state (DOS) and gives the model a wider 
applicability. 



6 Thermodynamic Properties 

The relation between then QLG and the anisotropic Heisen- 
berg model is such that the chemical potential [i corre- 
sponds to the external field h z , i.e. the grand canonical 
partition function in the QLG corresponds to the canoni- 
cal partition function in the anisotropic Heisenberg model 
where the number of spins is fixed. Consequently, thermo- 
dynamic potentials of the two models are related as: 



Oqlg -fJ,N = 6 



Heisenberg 



(22) 



Here refers to an arbitrary thermodynamic potential. 
In the same way as the ground state minimizes the in- 
ternal energy U = (H) at absolute zero, the free energy 
F — U — TS is minimized at finite temperatures. We wish 
to stress that the internal energy in the present approxi- 
mation cannot be derived accurately from the expectation 
value of the Hamiltonian in the following way: 



i ij 



(23) 



The cumulant decoupling, though a good approximation 
to the anisotropic Heisenberg model, is also an exact solu- 
tion of an unknown effective Hamiltonian H e f f . Therefore 
thermodynamically consistent results are only obtained if 
the effective Hamiltonian is substituted in the equation 
above, i.e. U = (H e ff). Here, as we do not know the ex- 
plicit form of this effective model, we have to integrate the 
free energy from thermodynamic relations: 



dF = 



(SI) + (S' B ) 



dh z + SdT 



(24) 



This equation allows us to select the ground state in re- 
gions where two or more phases exist self-consistently ac- 
cording to the random-phase approximation. The entropy 
of the spin system is given by 



S = J d 3 k 



1 



, 1 + exp(-/3w(k)) l0g( ^ 1 + exp(-/3w(k)) 



1 



+ 



1 + exp(/3w(k)) l0g<y 1 + exp(/3w(k)) 



) (25) 



This formula is of purely combinatorial origin and reflects 
the fact that the hard-core boson system is equivalent to a 
fermionic system given by the Jordan- Wigner transforma- 
tion. The w(k) terms refer to the energies of the spin-wave 
excitations. In the solid phases both branches have to be 
taken into account. The entropy of the anisotropic Heisen- 
berg model given for a fixed number of spins corresponds 
to the entropy of the QLG model at a constant volume. 
Therefore, in order to obtain the usual configurational en- 
tropy of the QLG, we have to divide by the number of 
particles per unit cell: (S conf = n J^ nB )- 
Here n A = 1/2 — S A and n B = 1/2 — S B are the particle 
occupation numbers on lattice sites A and B. As the na- 
ture of the normal solid to supersolid phase transition is 
not yet satisfactorily understood recent experiments [19] 
have focused on the behavior of the specific heat across 
the transition line in the hope of shedding light on the 
matter. The specific heat at constant temperature and 
constant pressure respectively are given by: 



Cx 



T 



dS, 



conf 



dT 



Cn = T 



T,V,N 



dS. 



'conf 



dT 



T,P,N 



(26) 

Although the external magnetic field h z in the spin model 
is an observable the corresponding quantity in the QLG 
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model, namely the chemical potential is not. Therefore we 
are also interested in attaining a formula for the pressure 
associated with a certain chemical potential. The relation- 
ship is most easily derived from the following Maxwell re- 
lation: 



dP 



T,V 



dN 
W 



# 



lattice— sites 



(l-e) (27) 



where e := (SV/ + (Sg). Note that, using this equation in 
order to obtain the pressure at any specific temperature 
we need a reference point, i.e. a chemical potential where 
the corresponding pressure is known. This point is given 
by /i — > oo which corresponds to N — ► and, hence, 
P — > 0. Consequently, in order to obtain the pressure for 
a specific chemical potential \i! we have to integrate over 
the interval [oo, //]. 



7 Boundary Lines 

The state of the system at any point in T and h z is given 
by the self-consistency equations and the free energy as 
can be derived from Eq. (f24|) . Nevertheless the resulting 
computations come at high computational cost and there- 
fore it seems most feasible to derive ODEs which deter- 
mine the first and second order transition lines. First we 
will derive the ordinary differential equations which define 
the more abundant second order transitions. The normal 
fluid (FE) and the normal solid (CFE) phases are deter- 
mined by Equations (l2Tj) and the supersolid (CAF) and 
superfluid (CFE) are defined by Equations (|20|) and con- 
dition Equation ©. Consequently on the SS-NS and SF- 
NF transition line, where both the normal (FE and CFE) 
and the super (CFE and CAF) phases coexists following 
equations hold: 



FZ((S tA ),{S, B ),h',T) = 
Fg((S ZA ),(S ZB ),h*,T) = 
(S* B ) 



h z + 2(S ZA )(4 - J 2 ' ) + 2{S ZB ) 4 = 2J X 



h z + 2{S ZB )(4 - Jj) + 2(S ZA )4 = 2J~l 



-(S ZA ) 



On the SS-NS boundary line the quotient (S XA )/(S ZA ) is 
not known a priori and therefore we eliminate it in the 
equation above yielding: 

FZ((S ZA ),(S ZB ),h* ) T) = 
F»((S ZA ),(S ZB ),h*,T) = 
f{(S ZA ),(S ZB ),h z ) 

:= (h' + 2(S ZA )(4(0) - Jj) + 2(S ZB )4) x 
(h* + 2(S ZB )(4 - J 2 T ) + 2(S ZA )4) - (2JJ{S ZB )) 2 =(20) 

We introduce a variable s which parametrizes the bound- 
ary curve. If we for instance choose ds=dT we get the 



following set of ordinary differential equations, defining 
the NF-SF and the NS-SS transition lines: 





0F% 


9{S, A ) 


d(Sz B ) 


9F B 


dF§ 


d(S, A ) 


9(S, B ) 


df 


Of 



9(Sz A ) 







OF™ dF% \ 
dh* dT • 

Oh* 
df 
Oh* 





dF± 
dT 



1 J 



ds 



dh ; 



\ ds I 



/cr 
o 
o 

Vi 



(30) 



Upon crossing the SF-SS transition line, coming from the 
superfluid phase the set of possible solutions branches off 
into two phases, the supersolid and a non-physical (com- 
plex valued) superfluid phase. Therefore any matrix of or- 
dinary differential equations will render a singularity and 
consequently we have to approach the transition line in a 
limiting process: 

Fi({S XA ),(S XB ),h*,T) = 

F§((S XA ),(S XB ),h*,T) 



lim((S XA ) - (S Xh 

e— >0 



The resulting ODE is 
/ 



9F% 
d{S ) 
dFg 



d{SZ 
1 



_9Fa 

d{S x 
OFi 
d^ 
-1 





dh' 
9F S E 
dh' 






dF± 
dT 

dF§ 
dT 



1 





e) = 

/ 9(S XA ) \ 
ds 1 



.9s 



V 



/ 



o 



(31) 



(32) 



which defines the superfluid to supersolid transition. As 
mentioned previously, there are certain parameter regimes 
where not all four possible phases are appearing and con- 
sequently a first order transition (mostly superfluid to su- 
persolid) will occur. In the previous chapter we have seen 
that such a transition line is difficult to locate. However, 
a tricritical point frequently appears in the phase diagram 
and at this point the first order transition evolves into a 
second order transition. This tricritical point can be taken 
as a initial value for a differential equation defining the 
corresponding first order transition line. 
The relevant ODE may be derived from a Clausius Clapey- 
ron type equation. On the transition line both phases have 
equal free energy. Hence: 

AS = K (33) 

A((S ZA ) + (S ZB )) dh z { 1 

where S refers to the spin entropy as derived in the pre- 
vious section (Eq. (|25p) and A refers to the difference of 
either entropy or spin mean-field of the superfluid and the 
normal solid phases. This equation together with ds = dT 
and the total derivative of the two self-consistency equa- 
tions for the normal solid and one equation for the super- 
fluid phase form a set of 5 ODEs determining (S x ) in the 
superfluid phase and (<5' ZA } and (S ZB ) in the normal solid 
phase along the boundary line in the T — h z plane: 

n §F± dF^ \ 

U V i^uz Pit 1 



1 dW* 








9FI 

d(S, A ) 

d{S z 






0*2 

d{S zl 
OFi 






gFf 

dh' 
d£l 
dh' 
d£l 
dh' 

AS 




dF b 
dT 
dFf 
dT 
d£l 
dT 

A(S ZA+B ) 
1 



(34) 
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8 Excitation spectrum 

The superfluid phase features, due to spontaneously bro- 
ken U(l) symmetry, the well know gapless Goldstone bosons, 
i.e. linear phonons. The supersolid phase additionally ex- 
hibits a second, gapped branch which is due to the break 
down of discrete translational symmetry. 
Figure 2] reveals a zero frequency mode in the supersolid 
phase at [100] of the first Brillouin zone and consequently 
the superfluid to supersolid phase transition is character- 
ized by a collapsing roton minimum at [100]. The disper- 
sion relation in the superfluid (CFE) phase is given by: 



c(fc) = 2{(J 1 T ( 7l (fc)- 
[<S z ) 2 (Ji T (7i(fc) 



1 



H (72(fc) 



1) + JlMk) 



l))x 
-1))- 



(S X ) 2 (JJ + J 2 T - Jhi(k) - 4l2{k))\ } 1/2 (35) 

From this equation we can see that the energy possibly 
goes to zero at [100] (corresponds to 71 (fc) = —1 and 
72 (&) = 1) when following condition is fulfilled: 



J{ + J2 



< 



(36) 



Hence we obtain a further condition (supplementary to 
Eq. (fT6|) ) for the existence of the superfluid to supersolid 
transition. 

Equation ([33)) allows for the existence of a second region 
of the reciprocal lattice space where the dispersion rela- 
tion might go soft. For 7i(fc) = and ^(k) = —1 which 
corresponds to [111] we obtain following condition: 



-2JJ 



J? 



> 



•Jo 



(37) 



It is also interesting to study the behavior of the excitation 
spectrum with increasing temperature. In a conventional 
superfluid the long wave-length behavior is given by: 



U (k) = n ° { ™ 0) k 



(38) 



Here V(0) is the interaction potential at zero momentum 
and m is the particles' mass. The density of the conden- 
sate n (T) typically decreases with increasing temperature 
and for that reason we expect lower energies with increas- 
ing temperature. In Figure [3] we can see that the quasi- 
particle energies indeed decrease with increasing temper- 
ature. Apart from the region in the vicinity of [100] the 
energies at higher temperatures lie significantly lower than 
those ones closer to absolute zero. This is important as it 
will contribute to the variation of thermodynamic quan- 
tities such as the entropy or the specific heat. Figure 0] 
depicts the variation of the excitation spectrum with in- 
creasing temperature in the supersolid phase. In this phase 
the low lying phonon branch mostly depletes with increas- 
ing temperature, although there exist a region between 
[110] and [111] where the zero temperature spectrum is 
significantly higher. Contrary to the phonon branch the 
gapped mode lifts the energy around the long wave length 
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[000] [100] [110] [111] [000] 

Fig. 3. Excitation spectrum in the superfluid phase for for 
Jj = 1.498if, Jj = 0.562K, jj 1 = -3.899K, 4 = -l-782Jf 
and h z = 3. Solid line refers to T = 0, dotted to T = 0.5, 
dashed line to T = 1 and the long dashed line to T = 1.3. 




[000] [100] [110] [111] [000] 

Fig. 4. Excitation spectrum in the supersolid phase for for 
Jj = 1.498A-, Jj = 0.562if, j] 1 = -3.899K, j\ = -1.782K 
and h z — 0.8. The Solid lines refer to T = and the dashed 
lines to T = 1. 



limit [000] and around [100]. In comparison with the super- 
fluid dispersion relation the excitation spectrum changes 
its form/shape rather than scaling down with increasing 
temperature as in the superfluid phase. We observe that 
the supersolid phase exhibits a more complex and diverse 
structure than the superfluid or normal solid phases alone. 



9 Discussion 

9.1 On Finite Temperature Properties 

In this section we will discuss finite temperature properties 
of the anisotropic Heisenberg model and its solution in 
the random-phase approximation. In order to be able to 
compare the temperature dependence of the model in the 
random-phase approximation with the classical mean-field 
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T[KJ 



Fig. 5. The phase diagram for J{ = J 2 ' = 0.5627f, 

jf = -3.899if and jj = -1.782-K' in the RPA (solid lines) 
and in classical M F (dashed lines) . MF overestimates the tem- 
perature by about 30%. 

approximation we chose the set of parameters that was 
extensively scrutinized by Liu and Fisher [12]: 

J7 = 1.4981s: 
J 2 T = 0.562X 
jf = -3.8991C 

4 = -1.7821s: (39) 

As mentioned in the previous section, Liu and Fisher [12j 
have chosen this set of parameters because it provides ar- 
guably the best fit to the phase diagram of Helium-4. Since 
we believe that the validity of the quantum lattice gas 
model is too limited to appropriately reproduce the be- 
havior of Helium-4 over the whole range of temperature 
and pressure we do not discuss most properties in the 
pressure-temperature space but rather present the major 
part of the results in the more comprehensible magnetic 
field -temperature coordinates. Only where the theory can 
be compared to relevant experimental data, such as the 
heat capacity at constant pressure we work in the corre- 
sponding representation. 

The phase diagram of the anisotropic Heisenberg model in 
T and h z coordinates is given in Figure[5l The dashed lines 
correspond to the phase diagram of the mean-field ap- 
proximation. We see that the diagrams are quantitatively 
quite similar but in the mean-field solution the tempera- 
ture is somewhat overestimated giving an approximately 
30% higher temperature for the tetra-critical point. As 
mentioned before the critical magnetic fields h z , due to 
quantum fluctuations, are lower in the random-phase ap- 
proximation. This effect is most distinctive on the super- 
solid to superfluid transition line as there the deletion of 
the spin magnitude is strongly pronounced. 
The net vacancy density e, density of vacancies minus den- 
sity of interstitials sparked interest as the question arose 
[TB] as to whether the number of vacancies follow predic- 
tions of thermal activation theory or are due to the effects 
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Fig. 6. Net vacancy density for four different pressures. The 
dashed line shows the curve expected if the normal solid would 
exist down to zero temperature 
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Fig. 7. Net vacancy density as a function of the external mag- 
netic field (equates to the chemical potential) for four different 
temperatures. The solid line refers to T=0K, the dashed line 
to T=0.4K, the long dashed line to T=0.7K and the dotted 
line refers to T=1K. 



of an incommensurate crystal. Figure [6] shows the net va- 
cancy density in the supersolid and the normal solid phase 
at constant pressures. In isobar curves, curves of constant 
pressure, the magnetic field h z is controlled through Equa- 
tion (|2T|) . We see that the net vacancy density is nearly 
constant in the supersolid phase. Only in the normal solid 
phase the net vacancy density increases exponentially with 
increasing temperature in total agreement with classical 
thermal activation theory. The almost temperature inde- 
pendent behavior of the net vacancy density in the su- 
persolid phase is an important finding of the quantum 
lattice gas model and should be observable in high reso- 
lution experiments if this effect is real. In Figure ([7|) we 
plotted the net vacancy density as a function of h z (or 
equivalently \i the chemical potential) across the normal 
solid, the supersolid and the superfluid phases for various 
temperatures, namely T = OK, T = OAK, T = 0.7K 
and T = 1.0K. The term net vacancy density does not 
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T/Tc 



Fig. 8. Free energy in the supersolid phase. The leading con- 
tribution comes from a T 4 ' 2 -term. The inset shows the free 
energy in double logarithmic scale. The dashed line is the fit 
to the T 4 ' 2 -term. The leading correction comes from a T 5 -term, 
indicated by the long dashed line. 




Fig. 9. [Color online] 3D plot of the spin entropy. The four 
phases are clearly distinguishable as the entropy is non-smooth 
across the transition lines. 
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Fig. 10. Specific heat at constant pressure. 
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e = s V s 'b 

Fig. 11. Curves show the net vacancy density as a function of 
temperature on the supersolid to normal solid transition lines. 
The solid line refers to set 1: jj = 1.498-ff, Jj = 0.562^, 
jf = -3.899A: and jj = -1.782.K' and the dashed line to set 
2: jj = 0.5K, Jj = 0.5K, jj 1 = -2.0K and jj = -0.5K. 

have a physical meaning in the superfiuid phase and here 
the quantity e rather refers to the particle density of the 
fluid. The particle density in the superfiuid phase is lin- 
ear in h z and independent of the temperature T, which 
follows immediately from condition Eq. ([9]) as (S%)/(Sg) 
is equal to one in the fluid phase. In the supersolid phase 
the dependence of the net vacancy density on the chemi- 
cal potential is stronger than in the superfiuid phase. The 
chemical potential (magnetic field h z ) is roughly inversely 
proportional to the pressure. Meaning that the superfiuid 
phase exhibits a higher compressibility than the supersolid 
phase, which is a quite remarkable result. Interestingly the 
net vacancy density in the supersolid phase also increases 
linearly with the chemical potential. This is due to the 
ratio of the superfiuid order parameter (S^)/(Sg) varies 
as the square root of the magnetic field in the vicinity of 
the transition: 

(£Z)/(S%) cx 1 - y/h> - h-CAF—AF (40) 

The exponent ^ is typical for mean-field type approxima- 
tions and appears close to all transition lines. Figure [TT1 
shows the net vacancy density of the model on the su- 
persolid to normal solid transition line as a function of 
temperature, and also reveals the mean-field type square 
root law dependence. 

At zero temperature in the solid phase the net vacancy 
density is equal to zero, hence the crystal is incompress- 
ible. In real systems compressibility usually occurs as a 
result of change in the lattice constant a, characterized by 
the Grueneisen parameter. The quantum lattice gas model 
does not take this effect into account since the lattice con- 
stant is treated as a constant. However measurements have 
shown that the lattice constant in the (supersolid) helium 
is almost a constant indicating that the net vacancy den- 
sity is the crucial parameter. Also of interest are the free 
energy and the entropy. Figure ([8|) shows the free energy of 
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the anisotropic Heisenberg model in the supersolid (CAF) 
phase at constant pressure. At low T the free energy usu- 
ally follows a power law: 



f oc r 



(41) 



The coefficient a is a universal property of the model 
which is constant over the entire regime of a phase. The 
exponent is most easily acquired in the log- log plot, shown 
in the inset of the plot. In this logarithmic scale a is given 
by the slope of the curve and the leading contribution is 
given by, approximately 



a = 4.2 



(42) 



in the supersolid region. This is close in value to the usual 
T 4 -term attributed to the linear phonon modes, and fol- 
lows from Equation (|25[) with u>(k) linear in k. Here the 
T 4 -term is solely due to the superfluid mode and in real 
solids an additional T 4 -term contribution, accounting for 
the lattice phonon-modes, will appear. The logarithmic 
plot also reveals the leading correction to the free energy 
given by a = 5. 

The non-configurational entropy of the system over the en- 
tire range of temperature and magnetic field h z is given in 
Figure O All four phases are visible and as expected from 
a thermodynamically equilibrated system the entropy is 
monotonically increasing with respect to temperature. 
Figure [TOl depicts the configurational specific heat at con- 
stant pressure in the supersolid and the normal solid phase. 
The jump in the specific heat at the critical temperature 
T c , in agreement with the second order phase transition, 
appears to be smeared out due to numerical inaccuracies 
as the specific heat is the second derivation of the free 
energy which had to be integrated of the interval [h z , oo). 
The jump in the specific heat may also be calculated from 
following formula which is an analogy to the Clausius- 
Clapeyron equation: 



A 



d 2 F fdh z 



dh* 2 V dT 2 



A 



TL 



d 2 F /dh z \ 
dh z dT VdTJ 



d 2 F 
A— = 



TL 



dT 2 



(43) 



AC h = -TA 
= TA 
+TA 



d 2 F 
dT 2 

d((S ZA ) + (S ZB )) 



dh z 

d((S ZA ) + (S ZB )) 



dT 



dT 

dhz 
dT 



TL 



TL 



(44) 



As we have Cp = Ch — T dT Q h (^t)p, we have for the 
specific heat at constant pressure: 



AC P = TA 
+TA 



d((S ZA ) + (S ZB )) 



dh z 

d((S ZA ) + (S ZB )) 



dK 
dT 



TL 



dT 



dT 



TL 



dh^ 
dT 




T [K] 



Fig. 12. The phase diagram for 3{ = 0.5K, jj = 0.5K, 
j] = — 1.0K and jj = 0.5K. The supersolid phase vanishes 
below T=0.323K. The resulting SF-NS transition is first order. 



For the values corresponding to Figure [10] we obtain an 
estimated jump of 0.02 which is in good agreement with 
the curve. 



9.2 First Order Boundary Lines 

In parameter regimes where the supersolid phase does not 
appear in certain temperature regions there consequently 
appears a first order phase transition between the super- 
fluid and the normal solid phase. Liu and Fisher [12j com- 
pared the free energies of the competing phases to estab- 
lish the transition line. This procedure is not applicable 
in the random-phase approximation, as was outlined in 
the section on thermodynamic properties. Other than the 
mean-field approximation where the Hamiltonian is given 
by Equation (J6j) the effective Hamiltonian of the random- 
phase approximation is not known. Therefore we have 
to integrate the first order transition line from a Clau- 
sius Clapeyron like equation as derived in the section on 
Boundary lines. A set of parameters which exhibits such 
a first order transition at low temperatures is given by: 

Jj = 0.5 
Jj = 0.5 



jf = -1.0 
4 = 0.5 



(46) 



(45) 



The corresponding phase diagram is shown in Figure 1121 
According to mean-field Eq. (fT4| and Eq. (|16p the sec- 
ond order superfluid to supersolid transition as well as 
the supersolid to normal solid transition is at absolute 
zero at h z = 0.86603, implying that the supersolid phase 
does not exist at zero temperature. At higher tempera- 
tures (above T > 0.323-ftT) the supersolid phase does exist. 
At T = 0.323X where the SF-SS and the SS-NS transi- 
tion lines intersect there occurs a tricritical point. On this 
tricritical point the superfluid and normal solid phases co- 
exist (as well as the supersolid phase) and the first order 
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tricrilical points 
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Fig. 13. The phase diagram for jj = 0.38K, jj = 2.5K, 
j\ = —1.2K and jj = 2AK. The supersolid phase does not 
appear above T=1.785K. The two tricritical points are con- 
nected by a first order SF-NS phase transition. 



4 



3 



2 









1 1 1 1 

NF 


1 


■ 


- 


SF 




- 












unknown phase 
1,1, 


, 





1 1 ' 1 ' 1 ' 1 ' 1 

0.2 0.4 0.6 0.8 

T[K] 



Fig. 14. Phase diagram for J X T = 0.5K,jJ = 0.5K, jf = 
—2K and jj = —1.5K. The superfiuid phase becomes unstable 
due to a imaginary quasi-particle spectrum (dashed line) . RPA 
does not yield any stable phase beneath that line. 



phase transition line can be calculated from the ordinary 
differential equation as has been derived in the section on 
boundary lines (Eq. {34])). The tricritical point is as such 
the starting point for the integration of the ODE. 
There also exist a regime of parameters where the super- 
solid phase is suppressed at higher temperatures as can be 
seen in Figure [13] which corresponds to: 

J X T = 0.38 

Jj = 2.5 
j\ = -1.2 

4 = 2.4 (47) 

The SF-SS and the SS-NS transition lines converge before 
the NF-SF line is reached, hence no tetracritical point 
such as in Figure [5] is present. The two tricritical points 
are connected by a first order superfiuid to normal solid 
transition line. 

9.3 Beyond the Model 

In Section [8] we have seen that the supersolid phase ap- 
pears if and only if the roton dip at 71 = — 1 and 72 = 1 , 
i.e. [001] of the first Brillouin zone collapses. Additionally 
we have seen that the model also allows for a collapsing 
minimum at [111] if condition Eq. ([37]) is met. A set of 
parameters that fulfills this condition is given by: 

Jj = 0.5K 
J 2 T = 0.5.K" 
jf =-2K 
4 =-1.5K (48) 

Note that the nearest neighbor and the next nearest neigh- 
bor constants in this configuration jj and jj , corre- 
sponding to the kinetic energy are relatively weak and 



are both of equal strength, leading to a highly anisotropic 
kinetic energy. Figure [H] shows the corresponding phase 
diagram according to the random-phase approximation. 
The normal fluid to superfiuid transition line starts at 

h z = 4.5 (49) 

at absolute zero and decreases with increasing tempera- 
ture. According to Eq. ([14]) and Eq. ([16]) the critical exter- 
nal fields h z defining the superfiuid to supersolid and the 
supersolid to normal solid transitions are, due to negative 
values under the square root, imaginary and hence phys- 
ically not relevant. Consequently in the classical mean- 
field approximation the superfiuid phase extends down to 
h z — and a first order superfiuid to normal solid tran- 
sition docs not occur as the relatively large negative jj 
increases the free energy of a possible solid phase. 
The random-phase approximation however draws a slightly 
different picture. Analogous to the classical mean-field 
solution the random-phase approximation also yields a 
phase transition near h z = 4.5. But unlike the classical 
mean-fields solution, the superfiuid phase here does not 
survive all the way down to h z = 0. Due to the partic- 
ular choice of parameters the superfiuid phase becomes 
unstable at around h z = 2; i.e the quasi-particle spectrum 
turns imaginary at 71 (If) = and 72 (fc) = —1 ([111]). In- 
terestingly beyond this line no other stable phase exists in 
the present approach; there is no set of spin fields (S XA ), 
(S XB ), (S ZA ) and (S ZB ) that solves the self-consistency 
equations or (|2Tj) . 

To understand why the present approach breaks down and 
how the phase below h z = 2 might look like we first inves- 
tigate the physical meaning of the collapsing roton mini- 
mum at [100] (equivalently [010] and [001]) leading to the 
supersolid phase consisting of two SC sub-lattices as al- 
ready analyzed in the previous sections: The roton dip at 
[100] corresponds to a density wave given by: 

cos(27ra;/a) + cos(27ry/a) + cos(27rz/a) 
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Fig. 15. Two dimensional projection of the lattice structure 
of a (super)-solid phase on three sub-lattices triggered by a 
collapse of the roton minimum at [111] of the first Brillouin 
zone. 

This density wave takes on value one on sub-lattice A and 
minus one on sub-lattice B, hence it reproduces the peri- 
odicity of the supersolid crystal. 

In the same way a collapsing roton dip on the main diag- 
onals [111] as given here refers to following density wave: 

cos(7ra;/a) cos(iry/a) cos(7rz/a) (51) 

This density wave yields zero on sub-lattice B and alter- 
natingly one and minus one on sub-lattice A (neighbors 
have opposite signs). Consequently this phase refers to a 
(super)-solid phase exhibiting three sub-lattices A, B and 
C where the mean-fields take on three different values (see 
Fig. IT5]> . It would be quite interesting to study the pos- 
sibilities and properties of such a phase and we leave as 
future work the extension of the present approach to ac- 
count for a three sub-lattice phase and the investigation 
its properties. 

10 Conclusion 

In this paper we extended the mean-field theory of the 
QLG model by Liu and Fisher [12] at finite temperatures 
and employed a random-phase approximation, where we 
derived Green's functions using the method of equation of 
motion. We applied the cumulant decoupling procedure 
to split up emergent third order Green's functions in the 
EoM. In comparison to the MF theory employed by Liu 
and Fisher [12] the RPA is a fully quantum mechanically 
solution of the QLG model and therefore takes quantum 
fluctuations as well as quasi-particle excitations into ac- 
count. The computational results show that these quasi- 
particle excitations are capable to render the superfluid 
phase unstable and thus evoke a phase transition. Quan- 
tum fluctuations account for additional vacancies and in- 
terstitials even at zero temperature and most interestingly 
the net vacancy density is altered in the supersolid phase. 
A potential shortcoming of the RPA is the lack of knowl- 
edge of the effective Hamiltonian and therefore the inter- 
nal and free energy. Usually the free energy is needed to 
compute first order transition lines as the ground state is 
given by the lowest energy state. We bypassed this obsta- 
cle by deriving a Clausius-Clapeyron like equation which 
defines the first order superfluid to normal solid transi- 
tion line. The entropy which is an input parameter of 



this equation is calculated from the spin wave excitation 
spectrum. The jump in the specific heat across the sec- 
ond order superfluid to supersolid transition line reveals 
information about the nature of the transition. However 
the specific heat is the second derivative of the free en- 
ergy and thus the jump is smeared out by the numerical 
calculations. Consequently we derived an equation which 
gives an optional estimation of the jump and is in good 
agreement with the numerical estimate. Most important 
our theory predicts a net vacancy density which, in the 
supersolid phase is significantly different from thermal ac- 
tivation theory. In the normal solid phase the net vacancy 
density roughly follows the predictions of thermal acti- 
vation theory, although quantum mechanical effects give 
a measurable contribution. Across the phase transition, 
however, in the supersolid phase the net vacancy density 
stays rather constant as T increases. 
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